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. We present a class of nonlinear Schrodinger equations (NLSEs) describing, in the mean field 

approximation, systems of interacting particles. This class of NLSEs is obtained generalizing 
£Sj ' expediently the approach proposed in Ref. [G.K. Phys. Rev. A 55, 941 (1997)], where a 

classical system obeying to an exclusion-inclusion principle is quantized using the Nelson stochastic 
quantization. The new class of NLSEs is obtained starting from the most general nonlinear 
classical kinetics compatible with a constant diffusion coefficient D — h/2m. Finally, in the case of 
s-stationary states, we propose a transformation which linearizes the NLSEs here proposed. 
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^ ' One among the most investigated topics in quantum mechanics regards the nonlinear generalizations of the 
Schrodinger equation. In the last decades many nonlinear extensions of the Schrodinger equation have been pro- 
posed in literature either to explore the fundamental aspects of quantum mechanics, with the usual linear theory 
Q [ representing only a limiting case, or to describe particular physical phenomenologics. 

Among the many nonlinear variants of Schrodinger equations proposed in literature, we recall the Bialynicki-Birula 
and Mycielski |l| in which appears the term —b In \ip\ 2 , the Guerra and Pusterla model [2j introducing the nonlinearity 
A|'0|/|'0| in order to save the superposition principle of the linear quantum mechanics, and finally, the Weinberg theory 

which proposes the introduction of homogeneous nonlinear terms in order to save, only partially, the superposition 
principle. 

t— I \ Many of the NLSEs proposed contain complex nonlinearities. For instance, the Doebner-Goldin equation 0,0 , be- 
' longing to the Weinberg class, were introduced as the most general NLSE compatible with the Fokker-Planck equation 
for the probability density. 

In addition, a large number of NLSEs have been introduced in order to study some phenomenologies in condensed 
\ matter physics. In Ref.s [f| 0} a NLSE with a nonlinearity of the type a\ \i(j\ 2 tp + a-i |^| 4 ip + i a 3 c^d^l 2 t/0 + {a-i + 

1 a 5 ) d x \tl>\ 2 ip, is introduced to describe a single mode wave propagation in a Kerr dielectric guide. Another example 
is found in Ref. || where the nonlinearity a\ \ij)\ 2 ip + i ci2 tp + i a-3 d xx ip + i 04 \tp\ 2 tp takes into account pumping and 
dumping effects of the nonlinear media and can be used to describe dynamical modes of plasma physics, hydrodynamics 
and solitons in optical fibers (Ref. and references therein). Finally, NLSEs are used to describe propagation of 
high power optical pulses in ultrashort soliton communication systems 0, incoherent solitons |Jjj> an d others. 

In Ref. |12j . starting from a classical system obeying to an exclusion-inclusion principle, a NLSE with complex 
nonlinearity was derived by using a direct generalization of the stochastic quantization method proposed by Nelson 
[l3j . The classical nonlinear kinetics imposed by the exclusion-inclusion principle was studied in Ref.s [l4l Hif . 
Successively in Ref. ^(| the problem of the most general nonlinear classical kinetics was investigated. This last 
one contains as a particular case the kinetics found in different physical systems, like among others, the Bosons and 
Fermions (generalized exclusion-inclusion principle), the quons and the system of particles obeying to the Haldane 
statistics jTJJ. The above nonlinear kinetics can be used to describe several physical phenomenologies like the Bose- 
Einstein condensation, the superfluidity and superconductivity and so on. 

In the present paper we discuss the problem concerning the quantization of a classical system obeying to the most 
general nonlinear kinetics compatible with a constant diffusional coefficient, by using the method of Ref. We 
recall that in Ref. 0] it has been presented another generalization of the stochastic quantization method which 
permits to obtain NLSEs containing only real nonlinearities. 

Let us consider a classical stochastic Marcoffian process in n— dimensional space 1Z, described by the following 
evolution equation: 



dp(t,x) = 
dt 



[ir(t,y -»x)-7r(t,x-fy)] (Ty , (1) 



where p(t, x) is the particle density function and 7r(f , x — > y) the transition probability from the site x to the site y 
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given by 14]: 

n(t, x -> y) = r(t, x,x - y) a(p(t, as)) /3(p(t, y)) . (2) 

The above definition for the transition probability implies a nonlinear kinetics and plays a fundamental role in the 
development of the theory exposed here. We observe that ir(t, x — > y) is given by a product of three factors. The first 
one, r(t,x,x — y) 7 is the transition rate and depends only on the starting x and arrival y sites during the particle 
transition x — > y. The second factor a(p(t, x)) is an arbitrary function of the particle population of the starting site. 
This function satisfies the condition a(0) — 0, because if the starting site is empty the transition probability is equal 
to zero. Finally the third factor 0(p(t, y)) is an arbitrary function of the arrival site particle population. For this 
function we have the condition (3(0) — 1 which requires that the transition probability does not depend on the arrival 
site if, in it, particles are absent. The expression of the function j3(p{t, x)) plays a very important role in the particle 
kinetics because can stimulate or inhibit the transition x — > y allowing in this way to take into account interactions 
originated from collective effects. 

The evolution equation (JTJ, taking into account the expression of the transition probability Eq. J5J), assumes the 
form: 



dp(t,x) 

at 



P(p{t, x )) J r{t, x + y,y) ot(p(t, x + y)) d n y 
n 

- a(p(t, x)) J r(t, x, y) 0{p(t, x - y)) d n y . (3) 
n 

Now we consider the particle kinetics in the first neighbor approximation and then we expand up to the second order 
the quantities r(t, x + y, y) a(p(t, x + y)) and (3(p(t, x — y)) in Taylor series of x + y and x — y respectively, in an 
interval around x, for y <c x. We obtain 

r(t, x + y,y) a(p(t, x + y)) = r(t, x, y) a(p(t, x)) + d. [r(t, x, y) a(p(t, x))] y i 

+ i 5, 0. [r(t, x, y) a(p{t, as))] y z y. + 0(y 3 ) , (4) 

f3(p(t, x-y))= P(p(t, x)) - d, (3(p(t, x)) y t + ^8 z 8 j f3(p(t, x)) y t V] + 0(y 3 ) , (5) 
being d i = d/dx i , with i = 1, • • • , n. 

After introducing the drift coefficient J(t,x) the diffusion coefficient D(t,x): 

J(t,x) = J yr(t, x, y)d n y , (6) 
n 

D(t, x) 8 i} = \Jy, V 3 Kt, x, y) d n y , (7) 
iz 

Eq. transforms to the following nonlinear partial differential equation [lfil Il7^ : 

dp 



at v 



(8) 



with p — p(t, x), a = a(p) and (3 = /3(p). Eq. (JSJ is a continuity equation for the density of particle d p/dt + 'V ■ j = 
in the configuration space, where the current results to be the sum of two terms: j — J dri f t + Jdiff- The first one 
3 drift — — {J + ^ D) j(p) with 7(p) = a(p) fl(p), is a nonlinear drift current. The second term j diS = —D(f3da/dp — 
a d (3/d p) V p is a nonlinear diffusional current different from the standard Fick current. In fact, the coefficient which 
multiplies V p is a function of p. Eq. |(SJ describes a class of nonlinear diffusional processes varying the functions 
a(p) and (3(p). 

In the following, we study a sub-class of the above processes for which the diffusion coefficient D is a constant and 
the diffusional current is given by the standard Fick form _7 diff = —D V p. Then the functional a and f3 obey the 
condition: 



(9) 
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From Eq. © we obtain: 



,{ P ) = [ 1 {p)T{p)] 1/2 , (3(p) = h( P )/T(p)} 1/2 



(10) 



where 



T(p) = exp 



dp' 
lip 1 ) 



(11) 



Varying j(p), Eq.s ©, (|10|l and describe the class of nonlinear kinetics which is the object of study of the present 
contribution. This class of kinetics includes also the standard linear one for which j(p) — p and then a(p) = p, 
P(p) = l. 

Following Nelson, we assume that the continuity equation ijSJ), which is a Fokker-Planck equation, describes a 
stochastic Markofhan process asymmetric in time 13]. For this reason we must consider both the forward (+) and 
the backward (— ) expressions of the particle current j = j(p) =p D V p where = - J (±) are the forward and 
backward velocities. The forward and backward Fokker-Planck equation assume the form: 

dt 



u (±) 7 ( p ) T D V p 







(12) 



We define now the current velocity as v = (it'+) + )) /2, assuming that can be derived by a potential v = V S/m. 
After summing the forward and backward Fokker-Planck equations, we obtain the following continuity equation: 



dp 
dt 



l{p) 







(13) 



In the stochastic quantization method the system is viewed as a classical Markoffian process where the particles are 
subjected to a Brownian diffusion with a diffusion coefficient given by D — h/2 m jlil l2f|. The mean acceleration for 
this classical stochastic process is defined as 



1 

a = — 

2 



dt 



(-) 



••-+<£) 



(+) 



u 



(-) 



(14) 



where the forward and backward derivative are defined us: 



(±) 



d_ 

dt 



u 



:±) • V±DA . 



(15) 



The dynamics is governed by the Newton law F = ma and supposing that the mean external force F is derived from a 
potential F — —V V, we can write VV = — m a. This last equation, if we take into account the expression of a given 
by Eq. (|14|1 and the expressions of the forward and backward velocities = VS/m±DVp/7(p), reproduces after 
integration the evolution equation for the potential S which is the following generalized Hamilton-Jacobi equation: 



as 

~dt 



2 m 



+ U q (p) + W(p) + V = 



where U q (p) = — (Ti 2 /2 m) A p 1 / 2 / p 1 / 2 is the quantum potential |U|23| and W, given by 



W(p) = fi(p)Ap + f 2 (p) (Vpf 



with 



flip) = T- 
4m 



1 



hip) = 



4m 



1 

P liP). 

1 djjp) 1 
7 (p) 2 Op 2 7 (p)^ 



1 

2p* 



(16) 

(17) 

(18) 
(19) 



is a real nonlinear quantity depending only on the density of particle p. 
We introduce now the complex function ip given by the ansatz: 



</>(*, x) = pit, x) 1 ' 2 



exp 



~ Sit, x) 



(20) 
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By taking into account Eq. (|20|l the particle current j — (V S/rn) j(p) can be written as: 



i ft 
2m 



lip) 



Vtp Vtp* 

tp tp* 



(21) 



The continuity equation l|13|) , the generalized Hamilton- Jacobi equation (|16fl together with the ansatz l|2U|l , allow to 
obtain the evolution equations for the function tp which is the following NLSE: 



i ft 



dtp 



■ A tp + V tp + A{p, j) tp 



dt 2m 

where A{p, j) = W(p) + iW{p, j) is a nonlinear complex quantity whose imaginary part is given by: 



2p 



P - lip) 
lip) 



and the real one takes the form: 



W(p) 



4 m 



Ap- 



1_ dj _L_ 

7 2 dp 2 7 2 



1 

27 



(22) 



(23) 



(24) 



We remark that Eq. (1221) defines a large class of NLSEs by varying the functional lip). 

In Ref. fl^ ]. with the purpose of introducing a generalized inclusion-exclusion principle which originates from the 
collective effects of a many body quantum system, it has been assumed lip) = p (1 + up). The corresponding NLSE 
obtained from Eq. 112211 becomes: 



i ft. 



dtp 
~dt 



ft ft 

Atp + Vtp - in — V 

2m 2p 



PJ 



4m 



1 + Kp 

2-np {Vpf 



tP 



Ap 

l + Kp^ 2p {l + npf 



tp 



where j, given by Eq. I|21|l . takes the form: 



~itp* Vtp- tpVtP*)(l + np) . 
2 m 



(25) 



(26) 



In Eq.s (|25|l and H2(jfl the parameter n G R takes into account of exclusive (k < 0) or inclusive (k > 0) effects due to 
the particle interactions of the system. 

Finally we consider the s-stationary states (j = 0) of the system described by Eq. 11221) for which W(/5, 0) = 0. 



The field tp for these states has the form tp = p 

EtP = 



1/2 



exp(— i Et/ft) and obeys the stationary NLSE: 



2m 



Aij> + Vip + W(p)ip 



(27) 



Eq. H27f) describes a C-integrable system which is linearizable by means of a transformation on the field tp |23|, |2J 
In fact, if we introduce now the new stationary wave function <p by means of the ansatz 



(p = cT{p) 1/2 exp(-iEtfft) , 
from Eq. I|27|) follows that the new field cp obey to the stationary linear Schrodinger equation 

ft 2 , 



Ei 



2 m 



(28) 



(29) 



The constant c in Eq. I|28|l is obtainable from the condition J tp* tp d n x — J <p* (pd n x — 1 so that the wave-functions 
tp and (p can be normalized simultaneously. The study of the s-stationary states of the nonlinear system obeying to 
Eq. H27J) is reduced to the study of the s-stationary states of the linear system (I29|l under the same external potential 
Vix) and with the same energy E. We remark that the linearization of Eq. (|27(l can be understood by taking into 
account that the wave-functions tp and (p satisfies the relation (T<p)/cp — (Ttp)/tp + M^d^l 2 ) where T = —{ft 2 /2 m) A 
is the quantum operator of kinetic energy. This last relationship means that the sum of the kinetic energy and of 
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the nonlinear term VFdV'l 2 ) is equal to the kinetic energy of the system described by the wave function (j>. This last 
system is linear because the nonlinear quantity VF(|'i/ , | 2 ) is hidden into this kinetic energy. 

In conclusion, in the present paper we have studied within stochastic quantization a system obeying to the most 
general nonlinear kinetics compatible with a constant diffusion. The evolution equation (|22[) obtained in this way 
describes a wide class of NLSEs with complex nonlincarities by varying the functional "f(p). Being j(p) an algebraic 
function of p the nonlincarities are derivative of the second order, with the real part depending on the density of 
particle p but not on the phases S. All equations belonging to this class conserve the total number of particle of 
the system, they are L r (l)-invariant and, in absence of the external potential V{x), they are invariant over space- 
time translation and rotation, but generally, are not Galilei invariant. Finally, in the case of s-stationary states, a 
transformation 7-dependent permits the linearization of this class of the NLSEs. 
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